Theory of Fulde-Ferrell-Larkin-Ovchinnikov state of superconductors with and 
without inversion symmetry: Hubbard model approach 
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We study Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state of superconductors with and without 
inversion symmetry based on the Hubbard model on the square lattice near half-filling, using the 
random phase approximation. We show that center of mass momentum Q tends to be parallel to 
x- or y-axis in the presence of inversion symmetry, while Q vector is likely to be perpendicular to 
the magnetic field in the absence of inversion symmetry. We also clarify that d + /-wave pairing 
is favored and the hetero spin triplet /-wave state is present in the FFLO state unlike state in the 
superconductors only with the Rashba type spin-orbit coupling (RSOC) originating from the broken 
inversion symmetry. The triplet /-wave state is enhanced by magnetic field and the RSOC. This 
stems from the reduction of the spin susceptibilities by the magnetic field and the RSOC. 



I. INTRODUCTION 

Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state was 
predicted more than fourty years ago^* 2 - and has been 
intensively studied. 3,4 FFLO state is characterized by 
the formation of Cooper pairing in a magnetic field 
with nonzero center of mass momentum Q: (k f , — k + 
Q I). Up to now, various properties of the FFLO 
state have been predicted^' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 Recently, 
a strong candidate of the FFLO state, CeCoIn5, has 
been discovered ^&±ld&l&2&2L2L2&21& This material 
meets necessary conditions for the existence of FFLO 
state, namely being in the clean limit j2£ two dimensional 
electric nature^ paramagnetically limited upper criti- 
cal field^ and d-wave superconductivu y 29 ' 30 ' 31 , and is 
also known as a strongly correlated system. However, 
all previous studies of FFLO state did not take into ac- 
count the electron-electron repulsion beyond Fermi liq- 
uid corrections j 32 i 33 ' 34 i 35 Moreover, most of the studies of 
pairing symmetry of FFLO state are based on the varia- 
tional method. 

Another aspect of the FFLO state which has not 
been so emphasized is that since magnetic field breaks 
time reversal symmetry, singlet and triplet pairings 
should be mixed in the FFLO state according to the 
Pauli's principle, as in a ferromagnet attached to a sin- 
glet superconductor (FFLO-like state is expected to ap- 
pear in the ferromagnet in ferromagnet/superconductor 
junctions) ^> 3 - r .i 3 ^> 3 ^^i^- Actually, it is predicted that 
FFLO state becomes more stable by mixing p-wave 

• • 42 

pairing^ 

Recently, it has been found that FFLO-like state 
can also appear in surface superconductivity or 
noncentrosymmetric superconductors i 43 ' 44 ' 45 ' 46 ! 47 
Starting from the discovery of heavy fermion 
superconductor CePtaSi^ the study of non- 
centrosymmetric superconductors has reccnlty 
become a hot topic in condensed matter 

physics 49-50.51. 52. 53. 54. 55. 56. 57. 58. 59. 60. 61. 62. 63. 64. 65. 66. 67. 68. 69. 70 

Due to the broken inversion symmetry, Rashba type 



spin-orbit coupling (RSOC) is induced in the non- 
centrosymmetric superconductors) 71 ' 72 and therefore 
spin-singlet pairing and spin-triplet pairing can be 
mixed in superconducting state. 73 However, the study 
of FFLO state in noncentrosymmetric superconductors 
is insufficient, especially from the viewpoint of strongly 
correlated systems. 

In this paper, we study FFLO state of superconductors 
with and without inversion symmetry based on the Hub- 
bard model on the square lattice near half-filling, using 
the random phase approximation (RPA). By solving the 
linearized Eliashberg's equations directly, we can elude 
the difficulty of variational method. We show that cen- 
ter of mass momentum Q tends to be parallel to x- or 
?/-axis in the presence of inversion symmetry, while the Q 
vector tends to be perpendicular to the magnetic field in 
the absence of inversion symmetry. We also clarify that 
d + /-wave pairing is favored and the hetero spin triplet 
/-wave state is present in the FFLO state in constrast to 
the superconductors only with RSOC. The triplet /-wave 
state is enhanced by magnetic field and the RSOC which 
stems from the reduction of the spin susceptibilities by 
the magnetic field and the RSOC. 

The organization of the present paper is as fol- 
lows. In section II, we explain our model and intro- 
duce Eliashberg's equations with the RPA. In section 
III, we present calculated results of the eigenvalue of 
Eliashberg's equations and the gap functions. In sec- 
tion IV, a summary of the results in the present paper is 
given. 



II. FORMULATION 



We consider the square lattice without inversion center 
in a magnetic field oriented to the x-axis. The Hubbard 
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model can be written as 

H = -y~)(2t (cos k x + cos k y ) + fi) c\ a c k(J 

k.a 

(h + A sin k y \ 

-Asinfc K ■ o- s , s /4 s c fes > + y^n fcT n fci .(l) 

J k 

Here, fc represents two dimensional vector. We set lattice 
constant to be unity. The first term is the dispersion 
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with Matsubara frequency w, 
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relation and the second term consists of the Zeeman term 
with the energy h and the RSOC with coupling constant 
A. The third term represents on-site electron-electron 
repulsion. 



Then, the bare Green's functions have the following 
form in the 2x2 spin space: 
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FIG. 1: The diagrammatic representation of the linearized 
Eliashberg's equations. Broken line and Oi(i = 1 — 8) rep- 
resent effective interaction and spin indices, respectively. In 
the RPA, we take into account the contributions from bubble 
and ladder types of diagrams, which impose some relations 
on spin indices as shown in this figure. 



The linearized Eliashberg's equations with RPA in the 
weak coupling approximation are described as (see also 
Fig. QJ 
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AA SS (fc) = — (G so - (g + Q, iw m ) G S(J , (-q, -iu> m ) (-T ss (k - q)) 

+G- sa (q + Q, iu> m ) G- sa > (-q, -iu> m ) U 2 xt a d S (k + Q))^™' (?) > ( 5 ) 
AA s ,_ s (fc) = -^ £ (G s ,(q + Q,iu; m )G- s ^(-q,-iu m )(-T 3 ^ s (k-q) + U\^ d (k + q)))A aa ,(q), (6) 

-r s , s (fc) = i[/ 2 Xc (*) + ^ 2 xs (fc) , (7) 

-r s ,_ s (fe) = -U+ \u 2 Xc (k) - \u 2 XS (k) - U 2 lm X s '- s (k) (8) 
I 

with s =t) I) center of mass momentum Q, and inverse X ^ — X 1 ^ an d X^ — (x^)* are satisfied, x^ an d X^ 
temperature (3. Here, xs and xc are spin and charge sus- 
ceptibilities at uj n = 0, respectively, which are obtained 
by Xs = X^ — Rex^ and xc — X^ + R- e X^- Note that 
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are given by 



X TT (fc) 




A = 1 - U 2 '- 2 



2C/Re X2 , 



(9) 
(10) 



and n denote the temperature, k-point meshes and the 
band filling, respectively. For example, the magnitude of 
A in CePt 3 Si is A ~ 0.3^ 



Center of mass momentum 
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xild and xJL are defined as 



\xL(k))-B{ \ 2 1 J< (13) 



B = (1 - C/xi) 2 - |t/X2 1 2 - Notice that 



TT 



' and 



Xlad 



(xladj are satisfied. 
In the RPA, we take into account the contributions 
from bubble and ladder types of diagrams^ In the above, 
and x stem from the bubble type of diagrams, while 

T T t I 

Xiad ano - Xiad originate from the ladder type of diagrams 
(see Fig. [5]). Within the RPA, self-energy corrections 
and the frequency dependence of effective interaction are 
ignored. However, we can grasp the essence of the physics 
by the RPA^S 
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FIG. 2: The diagrammatic representation of x"" an d Xlad ■ 
The summation over the internal spin indices, e.g., a" has to 
be taken. 

By solving the Eliashberg's equations, we can obtain 
the gap functions. We define singlet component of the 
pair potential and triplet one with S z — for a later 
convenience as 



A s = (A n -A iT )/2, 
A t = (A n +A iT )/2. 



(14) 
(15) 



FIG. 3: (color online) Eigenvalue A as a function of the center 
of mass momentum at A = 0. (a) h = 0.01. (b) h = 0.1. 
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FIG. 4: (color online) Magnitude of Q as a function of h. 
\Q\ — h and \Q\ = 2h are depicted for comparison. 



III. RESULTS 

In the following, we set t = 1 and parameters as T = 
0.01,7V = 256 x 256, U = 1.7 and, n = 0.9. Here, T, N, 



Let us first discuss the center of mass momentum 
Q = (Q x ,Qy) in unit of 7r/128 in order to study the 
stability of the FFLO state. The superconducting tran- 
sition temperature is determined by the condition A = 1 . 
Since A is a decreasing function of temperature, the state 
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FIG. 5: (color online) Eigenvalue A as a function of the center 
of mass momentum at h = 0.1 and A = 0.3. 
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FIG. 7: Eigenvalue A as a function of the center of mass 
momentum at Q x = 0, h = 0.1 and A = 0.3. 
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FIG. 8: (color online) Fermi surfaces at h = 0.1, with (a) 
A = and (b) A = 0.3. Arrows denote nesting vectors. 



FIG. 6: (color online) Eigenvalue A as a function of the center 
of mass momentum at h = and A = 0.3. 

with larger A is favored. Thus, we study A as a function 
of Q to obtain optimal Q. 

Figure [3] shows the eigenvalue A as a function of the 
center of mass momentum at A = 0. For small magnitude 
of the magnetic field, h — 0.01, the state with negligibly 
small Q is favorable as shown in Fig. [3] (a). For a large 



magnitude of the magnetic field, h = 0.1, the finite mo- 
mentum state with Q x — 8 and Q y = is favorable as 
shown in Fig. [3] (b). Therefore, the vector Q tends to 
be directed to x-axis. The magnitude of Q is given by 
|Q| ~ h. In order to study the validity of this relation, 
we show magnitude of the Q vector as a function of h in 
Fig. |4j The approximate relation \Q\ ~ h holds for wide 
range of h. Note that the same result can be obtained by 
exchanging Q x by Q y , Q x by -Q x , or Q y by -Q y due to 
the four fold symmetry of the square lattice. Thus, we 
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can conclude that the Q vector is likely to be parallel to 



x- or y-axisi&i&i^ 



Re A,, 



A similar plot in the presence of the RSOC is shown 
in Fig. where we choose h = 0.1 and A = 0.3. In this 
case, the vector Q tends to have a finite value (Q = (0, 6)) 
and be parallel to y-axis. The magnitude of Q is roughly 
given by \Q\ ~ h. Note that finete momentum state is 
unfavorable in the presence of the RSOC but the absence 
of the magnetic field as shown in Fig. [6j where we choose 
h = and A = 0.3. The state with Q and that with — Q 
are not degenerate due to the inversion asymmetry as 
shown in Fig. [Jj where A is plotted for all Q y . 

The direction of the vector Q can be understood by the 
structure of split Fermi surfaces. Figure [8] depicts Fermi 
surfaces at h = 0.1, with (a) A = and (b) A = 0.3. 
It is mentioned in Refs£i22 that the good nesting vector 
is important for the FFLO state as can be seen from 
the term G(q + Q)G(—q) in the linearized Eliashberg's 
equations. The summation over u> n in G(q + Q)G(—q) 
leads to the term 



/(6r(g + Q))-/(-M-g)) 
&(g + Q) + M-g) 



(16) 



with a, a' = ±. When this term has a large value, 
£,a{q + Q)£a'(— <l) > should be satisfied (otherwise the 
numerator of this term becomes negligible). Therefore, 
the nesting vector Q should connect nearest Fermi sur- 
faces and cannot be ~ (±7r, ±tt). Thus, we see that Q 
vector tends to be parallel to x- or y-axis for A = and 
to y-axis for nonzero A as described by arrows in Fig. [5] 
We can also understand the result for nonzero A as fol- 
lows. Since larger Fermi surface moves to the van Hove 
singularity (0, 7r) and the smaller one moves to the oppo- 
site direction with the increase of the RSOC as shown in 
Fig. Efb), the formation of Cooper pairing on the larger 
Fermi surface, which has a finite momentum parallel to 
y-axis, contributes dominantly to the superconductivity. 
Thus, Q vector is oriented to y-axis for nonzero A. With 
a similar analysis, we can expect that Q vector tends to 
be perpendicular to the magnetic field in the presence 
of the RSOC for arbitrary direction of magnetic field. 
For magnetic field, h — (h x , h y ), the dispersion relations 
become 



Then, 



£± = — 2t (cos k x + cos k y ) — \i 
±\/ (hy - A sin fe^ 2 + (A sin k y + h x ) 2 . 



l/Xh x | 



(17) 




is satisfied when split Fermi surfaces intersect with each 
other. This indicates that split Fermi surfaces move to 
the opposite directions which are perpendicular to the 
magnetic field. Hence, we see that Q vector tends to be 
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FIG. 9: (color) Real and imaginary parts of gap function A|| 
at h = 0.1. We take A = 0, Q x — 8 and Q y = in the upper 
figure and A = 0.3, Q x = and Q y = 6 in the lower figures. 
Solid lines represent Fermi surfaces. Yellow arrow represents 
typical scattering process. 



perpendicular to the magnetic field in the presence of the 
RSOC, as predicted in other works Md^^. 

It should be remarked that A is reduced as ft. or A 
increase. Hence, the superconducting transition temper- 
ature is reduced with the increase of them. This can 
be also predicted by the reduction of spin susceptibilities 
(see Fig. [TJand Fig. Eg) . 



B. Pairing symmetry 

Here, we study gap functions which are normalized by 
the maximum value of |A S | as a function of wave vector 
k. Note that = — and A s ,A t £ R are satis- 
fied due to the weak coupling approximation. It is well 
known that singlet c?-wave pairing is dominant over other 
pairings in the absence of magnetic field and RSOCi 76 i 77 

Figure [H displays real and imaginary parts of the gap 
function Af-f at h = 0.1. We take A = with Q x = 8 
and Q y — in the upper figure and A = 0.3, Q x = and 
Q y = 6 in the lower figures, where A takes its maximum 
value (see Fig. [3] and Fig. [5]). As shown in this figure, 
it has a /-wave symmetry. Note that the imaginary part 
of A || is absent for A = 0. We see that the magnitude of 



(18) A|| is enhanced by magnetic field and RSOC. 



We show the real part of triplet gap function A t in Fig. 
[ITJl We also take A — 0, Q x — 8 and Q y = in the upper 
figure and A = 0.3, Q x = and Q y — 6 in the lower 
figure. As can be seen from this figure, it has a /-wave 
like symmetry. Note that in the absence of magnetic field 
but the presence of the RSOC, A t vanishesj 52 i 66 We show 
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FIG. 10: (color) Real part of triplet gap function A t . We take 
A = 0, Q x = 8 and Q y — in the upper figure and A = 0.3, 
Q x = and Q y = 6 in the lower figure. Solid lines represent 
Fermi surfaces. Yellow arrow represents typical scattering 
process. 



the real part of singlet gap function A s in Fig. QT] We 
can find that it has a d-wave symmetry. 

We can understand the appearance of /-wave symme- 
try in triplet components by the structures of the spin 
and charge susceptibilities (xs and Xc)- They have peaks 
near (±7r, ±7t) as shown in Fig. [T2] According to Eq.JTJ, 
the gap functions tend to have the same sign during the 
scattering process. Similar discussion is also applicable 

to Xild and Xiad ( see Fi S s - 1131 and ED, and also Eq.©. 
Therefore, /-wave symmetry is favored. We show typi- 
cal scattering processes toward (tt, 7r) by yellow arrows in 
Fig. 1 and Fig. [TO] 

Let us explain the origin of the enhancement of the 
triplet pairing. We study xs, Xc, x]l d and xJL^ which 
are plotted in Figs. E2rEU 6 - As shown in Fig. El xs is 
reduced by magnetic field and RSOC while xc is almost 
independent of them. This can be intuitively interpreted 
as follows. Magnetic field and spin-orbit coupling cause 
spin flip process and hence break magnetic fluctuation. 
On the other hand, spin flip scattering does not affect 
charge fluctuation. Therefore, xs depends on magnetic 
field and the RSOC while xc is almost independent of 
them. Since the positions of the peaks in xs and Xc 



ReA v 




FIG. 11: (color) Real part of singlet gap function A s . A = 0, 
Qx = 8 and Q y = in the upper figure and A = 0.3, Qx = 
and Q y — 6 in the lower figure. Solid lines represent Fermi 
surfaces. 



are almost the same, they compete with each other (see 
Eq.®). The decrease of xs leads to the reduction of 
singlet d-wave pairing and hence triplet pairings could 
dominate.— As shown in Fig. I13| Xiad' which contributes 
to the effective interaction for the heterospin pairings, is 
also reduced by magnetic field and the RSOC by the same 
reason, while the real part of Xiad> which contributes to 
the effective interaction for the equal spin triplet pairings, 
is enhanced by magnetic field and the RSOC. This results 
in the enhancement of the triplet components. Note that 
the imaginary part of x\ad emerges at h = 0.1 and A = 
0.3, although it is negligibly small except for the case in 
the presence of both of magnetic field and the RSOC. 

Next, we discuss the spin susceptibility for the re- 
direction, xT = \bdld + Re xJad)> and that for the v- 

direction, x V s = hi^ild ~ ^ e xjad)- The susceptibility for 
the z-direction Xs is defined as x| z = \Xs- We see that 
magnetic field strongly suppresses Xs x , while the RSOC 
strongly suppresses xf and xF : X V S V ~ Xs > XT at 
h = 0.1 and A = 0, and x| x > Xs' ~ Xs* ^ h — 0.1 and 
A = 0.3, as can be seen from Fig. [T21 and Fig. [15j Since 
the spin susceptibilities have peaks near (±7r, ±7r), the 
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FIG. 12: (color online) Spin and charge susceptibilities (xs 
and Xc)- We set h — 0.01 and A = in the upper figures, 
h — 0.1 and A = in the middle figures, and h = 0.1 and 
A = 0.3 in the lower figures. 



FFLO state coexists with antiferromagnetic fluctuation. 



TABLE I: Pairing symmetry and center of mass momentum. 
Here, n = (0, 0, 1), a; and y are unit vectors oriented to x- and 
y-direction, respectively, and h is magnetic field. 





h = A = 


h = 0,\=£0 


h=£0,\ = 






=0 


/-wave 


/-wave 


/-wave 


At 


=0 


=0 


/-wave 


/-wave 


A s 


d-wave 


d-wave 


d-wave 


d-wave 


Q 


=0 


=0 


II x,y 


| n x h 



Finally, let us summarize the main results in Table 
HI including results with h = A = and h = 0, A ^ 
0j 66 i 76 ' 77 The appearance of d + /-wave pairing in FFLO 
state is specific to strongly correlated electron systems 
(cf Refifi). 



IV. CONCLUSIONS 

In this paper, we studied Fuldc-Fcrrcll-Larkin- 
Ovchinnikov state of superconductors with and without 
inversion symmetry, where we used the Hubbard model 
on the square lattice near half-filling with the RPA. We 
clarified the following points: 




FIG. 13: (color online) xlld an d rea l part of xjad- We set 
h = 0.01 and A = in the upper figures, h — 0.1 and A = 
in the middle figures, and h — 0.1 and A = 0.3 in the lower 
figures. 




FIG. 14: (color online) Imaginary part of Xitd — 0.1 and 
A = 0.3. 



1. center of mass momentum Q tends to be parallel to 
x- or y-axis in the presence of inversion symmetry, while 
Q tends to be perpendicular to the magnetic field in the 
absence of inversion symmetry. This can be understood 
by nesting vectors of the Fermi surfaces. 

2. d+ /-wave pairing is favored. The hetero spin triplet 
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/-wave state is present in the FFLO state unlike state in 
the superconductors only with the RSOC. The triplet /- 
wave state is enhanced by magnetic field and the RSOC. 
This stems from the reduction of the spin susceptibilities 
by the magnetic field and the RSOC. 
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